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Abstract 

Atomic electrons are sensitive to the properties of the nucleus they are bound to, such as nuclear mass, 
charge distribution, spin, magnetization distribution, or even excited level scheme. These nuclear parame- 
ters are reflected in the atomic transition energies. A very precise determination of atomic spectra may thus 
reveal information about the nucleus, otherwise hardly accessible via nuclear physics experiments. This 
work reviews theoretical and experimental aspects of the nuclear effects that can be identified in atomic 
structure data. An introduction to the theory of isotope shifts and hyperfine splitting of atomic spectra is 
given, together with an overview of the typical experimental techniques used in high-precision atomic spec- 
troscopy. More exotic effects at the borderline between atomic and nuclear physics, such as parity violation 
in atomic transitions due to the weak interaction, or nuclear polarization and nuclear excitation by electron 
capture, are also addressed. 
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I. INTRODUCTION 



Atomic and nuclear physics historically share a common ancestry. When Niels Bohr proposed 
in 1913 his first model of the atom, he depicted it as having a small and dense positively charged 
nucleus, surrounded by the orbiting electrons. The existence of the electron, a discrete unit of 
negative charge, had been discovered in the cathode rays experiments performed by J. J. Thomson 
in 1897. However, it was the great interest arising from the phenomenon of radioactivity that 
engaged scientists in the quest for the microscopic world of the atom. The beginning of the XXth 
century was a very exciting time for physics, with one experimental discovery following the other. 
Quantum mechanics, the novel theory appropriate to describe the microscopic world of the atom, 
emerged to a good extent from the need to understand atomic spectra. With time, atomic physics 
turned to be the field primarily concerned with the arrangement of electrons around the nucleus 
and the processes by which these change, i.e., the atomic transitions. Nuclear physics, on the 
other hand, focused on the atomic nucleus, its building blocks and the interactions between nuclei. 
The question what does the nucleus consist of has led over the years to the Standard Model, the 
integrated theory of elementary particles, which is supposed to explain all processes in Nature 
ruled by the first three of the four known fundamental forces: the electromagnetic, weak, strong 
and gravitational forces. 

How much information about the nucleus is used in or can be obtained from atomic physics? 
This work intends to provide an up to date answer to these questions. Traditionally, atomic physics 
could afford to consider the nucleus as a point-like charge without loosing in accuracy. With in- 
creased precision in atomic physics experiments, the nucleus is going through an odyssey from the 
simplified picture of a point-like charge to its real size and properties. Different nuclear masses, 
charge distributions or spins do have an effect on the electronic levels and transition energies, 
and the small energy corrections or splittings observed in experiments can give in turn informa- 
tion about the atomic nucleus. Atomic spectroscopy, at present at an unrivalled level of precision 
achieved by frequency-comb techniques, can thus be used to determine nuclear parameters which 
are otherwise hardly accessible via nuclear physics experiments. The easiest way to determine 
the nuclear effect related differences in the atomic structure is to compare spectra of different 
isotopes — nuclei with the same atomic number, i.e., charge, but different mass number — and de- 
termine the small variation of electronic transition energies, the so-called isotope shifts. These 
shifts can be caused only by nuclear properties other than the charge. 
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The correct identification of nuclear effects in atomic transitions is very important in connec- 
tion with accurate theoretical atomic structure calculations involving many-body correlations or 
quantum electrodynamics (QED) corrections, which can be tested in high-precision experiments. 
A key ingredient for the determination of nuclear effects in atomic transitions and in general for 
tests of theoretical atomic and nuclear models are highly charged ions, advantageous due to their 
simplified electronic structure. Experiments at facilities with stored or trapped highly charged ions 
have played in the last decades an important role in the study of nuclear effects and determination 
of nuclear properties from atomic structure data. 

The borderline between atomic and nuclear physics is also inhabited by several nuclear pro- 
cesses that directly involve atomic electrons. The underlying interactions for such processes are 
either electromagnetic or weak. For instance, the nucleus and the electron can interact through 
the electromagnetic field and undergo transitions simultaneously. Such a process is well-known 
in nuclear physics, as being sometimes the only decay channel for a nuclear excited state: the 
internal conversion. A nucleus in an excited state that for some reason cannot decay radiatively 
transfers its energy through the electromagnetic field to one of the atomic electrons which leaves 
the atom. The inverse process of internal conversion, which might occur in highly charged ions, is 
called nuclear excitation by electron capture. The weak interaction, on the other hand, can create 
or annihilate atomic electrons, or can mix the parity of atomic states, leading to parity-violating 
effects. These processes at the interface between atomic and nuclear physics also belong to the 
category of nuclear effects in atomic transitions. 

In the following, theoretical and experimental aspects of these effects are addressed, structured 
in three parts. The first section deals with isotopic shifts, their origin in the charge volume and 
mass properties of the nucleus, and some of the established methods to determine them experi- 
mentally. The second section addresses the hyperfine structure of atomic spectra, its relation with 
the determining nuclear parameters, such as the nuclear spin and the nuclear multipole moments, 
as well as the various experimental methods used for observation. A separate section is dedicated 
to more exotic aspects such as parity violation and nuclear processes involving atomic electrons. 
The paper concludes with a discussion of the future directions in probing nuclear effects in atomic 
structure. 
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H. ISOTOPE SHIFTS 



By far the most important nuclear property that determines the atomic energy levels is the 
nuclear charge. The energy of the atomic levels can be obtained in the first approximation by 
solving the Schrodinger equation for the electrons in the central potential of the nuclear charge 
+Ze, where Z is the atomic number. In order to take into account the spin of the electron, the 
Dirac equation is introduced. However, in the first approximation, both equations describe the 
nucleus as a point-like charge of infinite mass. 

Measurements of atomic spectra for different isotopes of the same atomic number Z show slight 
differences — the isotope shifts. This frequency difference of an electronic transition is usually 
described separately as due to the finite mass of the nucleus — the mass shift — and to the size of 
the nuclear charge distribution — the volume of field shift. The mass shift is dominating isotope 
shifts for light atoms, whereas the field shift scales as Z 5 or even up to and thus becomes the 
leading contribution in the case of heavy ones. In order to extract information about the nuclear 
properties, a proper interpretation of the measured isotope shifts based on theoretical calculations 
is necessary. 

Since often nuclei with the same Z but different mass number A also have different nuclear 
spins, spectra of different isotopes may also have a completely different hyperfine structure. We 
will focus in this section on the mass and field isotope shifts, which do not give rise to splittings 
of the energy levels, and address the hyperfine structure in the next section. An introduction to 
the theory of the mass and field shifts is given in separate subsections, followed by an overview of 
several experimental methods used for isotope shift measurements. 

A. Mass shift 

Originally assumed for the sake of simplicity to be infinitely heavy, the nucleus obviously 
has a finite mass. The approximation of an infinitely heavy nucleus is acceptable for certain 
cases, but for a precise description one should consider the effects of nuclear motion or recoil. 
Nonrelativistically, in the case of hydrogen-like ions, this can be done easily by introducing the 
reduced mass, (see, for instance, Ref. yj]), 



where m is the mass of the electron and M the one of the nucleus. The Schrodinger equation is 
then solved in the center-of-mass system, with the reduced mass of the electron with respect to 
the nucleus replacing the electronic mass m. For many-electron systems, apart from the reduced 
mass, also a mass polarization term enters the Schrodinger equation, of the form 
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where the indices i and j denote the ith and jth electrons, respectively. Hughes and Eckart |2|] in 
1930 were the ones to separate the mass shift of atomic spectra into these two parts: the normal 
mass shift (NMS), due to the introduction of the reduced mass in the kinetic energy term in the 
Schrodinger equation (also called the one-body part), and the specific mass shift (SMS), which 
rises from the mass polarization term, or two-body part, introduced in Eq. ©. For decades, ac- 
cording to this separation, the NMS and SMS have been treated differently in calculations - the 
NMS was evaluated exactly, while the SMS was obtained from perturbation theory. This method 
of evaluation proved to be unreliable even in the non-relativistic case, as pointed out in Ref. yfl. 
Also, rigorously speaking, the mass shift theory in the nonrelativistic case is valid only for a point- 
like nucleus [|4|]. 

In the relativistic theory, NMS and SMS both rise naturally from the relativistic nuclear recoil 
operator, which can be written in the lowest-order relativistic approximation as 
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where r*j and pi are the coordinate and the momentum of the ith electron, respectively, di are 
the vector Dirac matrices, acting on the electron i, and a the fine structure constant. The NMS 
correction is given by the expectation value of the Rij operator for i = j, (J^ Ru), while the 
SMS correction can be obtained from the i ^ j terms, ^i)> treating the two terms on equal 

footage [3]. The full relativistic theory of the nuclear recoil effect can only be formulated in the 
framework of QED, see [la tZD- A precise optical isotope shift measurement with stored highly 
charged argon [|8|] ions has confirmed experimentally the necessity to consider the relativistic 
nuclear recoil effect in theory. More details about the experimental approach are given in Section 
III CI where the typical experimental methods to determine isotope shifts are outlined. 
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B. Volume or field shift 



Although often treated as point-like, the nucleus has a finite size. Consequently, the distribution 
of the protons in the nuclear matter determine the electrostatic nuclear potential, which will no 
longer be 1/r. The Dirac equation should be then solved numerically using a realistic nuclear 
potential instead of the usual Coulomb potential generated by a point-like charge +Ze. Compared 
to a point-like nucleus, the extended nuclear charge distribution leads to a shift of the resonance 
energies on the order of 200 eV in the case of the if-shell electrons in uranium (Z = 92). 

Since different isotopes of the same element with nuclear charge +Ze have different numbers 
of neutrons, their charge distribution will not coincide, leading to a small shift in their electronic 
level energies - the volume or field shift. Coming back to our previous example, in the case of the 
K-shell electrons in uranium, the field shift between 235 U and 238 U is less than 1 eV. The effect 
of the finite nuclear size is strongest for the s electronic orbitals, which have the largest overlap 
with the nucleus, and among these the field shift for the Is electrons is obviously the largest. The 
nuclear size effects scale with the nuclear charge as Z 5 up to Z 6 [9], such that the field shift value 
for the K-shell of uranium should be one of the largest measured field shifts. Uranium (Z=92) 
and californium (Z=98) are among the heaviest nuclei used for atomic physics experiments, both 
of them long-lived radioactive nuclei (in atomic physics, experiments with stable or long-lived 
isotopes are preferred due to technical difficulties). 

The history of isotope shifts goes back to 1931, when experiments by Schiiler and Keyston 
on the hyperfine structure of thallium and later on mercury Jiol U] led to the discovery of a 



structure which was not due to the nuclear spin but to a displacement of the atomic levels in 
different isotopes. The authors pointed out that their observations have to be explained by some 
differences in the nuclear fields of the isotopes. Ever since, several models for the nuclear charge 
distribution have been used in the literature. One of the obvious choices is to consider the nucleus 
as a homogeneously charged sphere of radius R. In this model, the electrostatic potential V(r) due 
to the nucleus is then given by yj] 

V(r) ={ 2R \ R 2 J (4) 
[ r>R. 

The nuclear charge distribution is in this case constant, p(r) = 3Ze/ (AnR 3 ). Using the nuclear 
charge distribution, radial moments (r n ) can be defined (see, for instance, Ref. Ill2ll ). of which the 
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most important is the root mean square charge (RMS) radius of the nucleus, 



<r 2 > 
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It is actually the RMS radius (r 2 ) 1 / 2 that the nuclear size effect depends on, rather than the charge 



distribution itself. Indeed, isotope shift measurements are (for nowadays precision) basical 



sensitive to the nuclear charge distribution, as long as the same RMS radius is reproduced 112 
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That is why most experiments aiming at the determination of the nuclear charge distribution extract 
the RMS radius out of the measured values with the help of various models. For the homogeneous 
sphere nucleus model, the RMS radius is given by (r 2 ) = 3_R 2 /5. 

A more physical description of the nuclear charge distribution is given by a two-parameter 
Fermi-type distribution |15i1 

p( r ) = ; ; P H_a/„ > ( 6 ) 



with the surface thickness 
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(7) 



The quantity c is the half-density radius which gives the distance to the point at which the nuclear 
charge density is one-half of its maximum value. The distance in the surface layer over which 
the nuclear density decreases from 90% to 10% of its maximum is called surface thickness and 
this is denoted by t in the above equation. A plot of the Fermi distribution highlighting the two 
parameters is presented in Figure [Q The Fermi distribution with t = reduces to a uniform 
distribution with radius R = c. A typically used value for all nuclei is t = 2.30 fm. The constant 
Po can be obtained by writing the total nuclear charge Ze as the integral over the Fermi-like charge 
distribution. Alternatively, the Fermi distribution parameters c and p can be expressed with the 
help of the RMS radius [14.H, 
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A commonly used shortcut to the nuclear charge distribution problem is that in most calculations 
the thickness parameter is taken as t = 2.30 fm, and the RMS radius is obtained from the empirical 
formula II 1511 

(r 2 ) 1 / 2 = (0.836A 1/3 + 0.570) fm , (9) 
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which is a weighted fit to the measured RMS radii for nuclei with A > 9. Its accuracy was found 
to be better than 0.05 fm and is in most cases sufficient for the present experimental accuracy. 

The two models for the nuclear charge distribution described so far assume a spherical nucleus. 
It is however well known that most of the isotopes have a deformed nuclear ground state. To 



account for the nuclear deformation, a deformed Fermi distribution can be introduced Iil7ll 



p(r,^) = 1 + e P ( ;_ e)/a , (10) 

where this time we have 

too I \ 
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In the equation above, f3i m are the nuclear multipole deformation parameters and Y im the spheri- 
cal harmonics. In most cases, axial symmetry is assumed, and considering only quadrupole and 
hexadecapole nuclear deformation, this expression reduces to 

c = c[l + (3 20 Y 20 {6, <p) + /3 40 Y 40 (9, tp)) . (12) 

This deformed nuclear charge distribution has been taken into account to calculate field shifts 
and derive approximate analytical formulas for nuclear- size corrections to binding energies of 



uranium and neodymium isotopes | |18|1 . The authors of Ref. [18] conclude that, if a precise value 
of the RMS radius is available from experiments, for an accurate calculation of the finite nuclear 
size effect on the electronic binding energies the nuclear deformation should be taken into account. 
The value of the finite-nuclear-size correction obtained from their analytical approximative method 
for H-like uranium 238 U 91+ is 198.39 eV when taking into account the nuclear deformation and 
198.61 eV when using a spherical nucleus with the same RMS radius. The theory thus predicts a 
220 meV energy shift between the spherical and deformed nucleus models. Whether it is possible 
to discern between nuclear models on the basis of such a small difference relying on experimental 
results is a different matter. In the following we address this practical issue and present some of 
the experimental methods used to determine the isotope shifts. 



C. Experimental determination of isotope shifts 

Isotope shifts are typically very small, and can only be determined in high-precision experi- 
ments. The better the experimental precision, the more one can learn about the nucleus and about 
the theoretical formalism that we use to obtain electronic transition energies. Usually atomic 
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transition energies do not directly provide information about the nucleus (other than the atomic 
number Z). Other parameters have to be extracted and interpreted by first taking into account all 
possible contributions to the electronic transition energy, whether from electron correlation, first- 
and second-order QED effects, nuclear recoil, nuclear finite size or nuclear polarization. Measured 
isotope shifts are evaluated in theory as the sum of the mass and field shift, plus QED and nuclear 
polarization contributions. Among these often the nuclear polarization contribution (which is the 
subject of Section HVBI) is the one with the largest theoretical error bars. As a bottom line, if 
one wants to extract out of the measured isotope shifts the most interesting nuclear property, the 
difference in the RMS radii of the different isotopes, one has to make sure that (i) the field shift is 
obtained out of the isotope shift by taking the correct theoretical values for all other contributions 
and (ii) a model for the field shift calculation delivers the RMS radius values. Thus, a reliable 
interpretation of the data is just as important as the measurement of the isotope shift itself. 

At the present level of accuracy, atomic physics measurements provide important information 
about the nucleus in the form of the already introduced RMS radius. The absolute values of the 
RMS nuclear radius can be obtained by two experimental methods: electron scattering and muonic 
atom x-rays. In both cases the theoretical analysis of the results is crucial for obtaining the nuclear 
parameter. 

Direct isotope shift measurements only provide the difference between RMS radii of two iso- 
topes, 5(r 2 ). Usually either the K a x-ray isotope shifts, or optical isotope shifts are used to deter- 
mine the variation of the RMS radii. Optical isotope shift measurements are performed by laser 
spectroscopy, a highly- sensitive method involving laser excitation of the transition of interest and 
subsequent detection of the decay. The difficulty here consists of interpreting the measured isotope 
shifts, especially due to the fact that the method is limited to low-charge states, i.e., many-electron 
systems. Atomic many-body effects make the theoretical description difficult and bring large un- 
certainties due to SMS and QED contributions in the calculation of the electronic density at the 
nucleus. X-ray spectroscopy can be used for the K a isotope shifts in highly charged ions. Less 
electrons make the theoretical interpretation easier, with the disadvantage that the experimental 
precision (one needs meV precision for x-ray transition energies) is limited. Alternatively, dielec- 
tronic recombination measurements on highly charged ions for a number of nuclei can achieve 
nowadays the necessary precision and provide complementary data to isotope shifts and RMS 
radii obtained via other methods. A brief description of the five methods, with special focus on the 
new and less traditional dielectronic recombination measurements will be given. Apart of dielec- 
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tronic recombination isotope shift measurements, which are relatively new, the RMS radii obtained 
from the other four methods have been combined in an attempt to elude the model-dependencies 
in the evaluation of the data in several compilations, from which the most recent is Ref. [19D. For 
an insight in the combined theoretical analysis of the four methods, see also II 1211 . 



1. Elastic electron scattering 

Elastic electron scattering refers to electron scattering in the Coulomb field of nuclei. Electrons 
are accelerated up to several hundreds MeV in a linear accelerator and then focused on a nuclear 
target. This energy is required in order to achieve a small wavelength for the electron despite its 
small mass. The electron beam should also be very intense, because of the very small scattering 
cross section, and have a very good energy resolution (unfortunately, these two criteria are not eas- 
ily achieved simultaneously, leading to the use of the so-called energy-loss spectrometer systems 



112011 ). What is measured is the differential cross section da(E } 9)/dVt for the elastic scattering of 
an electron of energy E at the angle 6. The angular differential scattering cross section is then 
compared with the relativistic formula of Mott for elastic scattering of electrons (see, for instance, 



a nice introduction of the scattering cross sections and form factors in Ref. [ 21(1 ). The Mott formula 
(in relativistic units for which A-keq = 1), 

da _ Z 2 e 4 cos 2 6/2 

dtt ~ 4£ 2 sin 4 fl/2 ' ( } 

is the relativistic equivalent of the well-known Rutherford formula for scattering of charged point- 
like particles on a charged point-like center. The finite size of the nucleus is introduced via a 
so-called form factor. For a spin-zero nucleus, in the first Born approximation, the scattering cross 
section can be factorized in (da / dVL) Mott for a point-like nucleus, and the form factor F(q), a struc- 
ture function which depends on the momentum transfer in the scattering, q = (2E/hc) sin(6 l /2), 
provided that m e c 2 <C E, 

\F(q)\\ (14) 
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The form factor is related to the nuclear charge distribution p„(r) through a Fourier-Bessel trans- 
formation tl2D 

F(q) = J d 3 rp n (r)j (qr), < q < oo , (15) 

where jo(qr) = (sin qr) / (qr) represents the spherical Bessel function of zero order. This Fourier- 
Bessel analysis has been introduced by Dreher et al. 12211 and provides a model-independent 
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method to obtain information about the nuclear charge distribution in elastic electron scattering 
experiments and in particular to extract the RMS radii. The first Born approximation is only suit- 
able for scattering on light nuclei, while for medium and heavy nuclei a careful analysis of the 
experimental cross sections should be done taking into account the phase shift J2I]. A number of 
RMS radii compilations for stable isotopes from electron scattering experiments and sometimes in 



19, 



21, 



2j. 



conjunction with other methods are available in Refs. H12I . 

In the last half a century, following the pioneering experiments by Hofstadter et al. Il21ll at 
Stanford, USA, electron-nucleus scattering has proved to be an excellent tool for studying the 
nuclei. Among the advantages of the method we have the generating electromagnetic interaction, 
which is on the one hand well understood, and on the other hand weak enough for the electrons 
neither to perturb the nucleus, nor to get absorbed or undergo multiple scattering. Also, this is 
the only electromagnetic method which can provide information about the radial dependence of 
the nuclear charge density p n (r). It should be noted, however, that the accuracy of the method 
is limited, with typical accuracies of 0.2% for RMS radii. Also, typical experiments work only 
with stable nuclear targets. With the development of Radioactive Ion Beam facilities, properties 
of exotic short-lived nuclei can become object of study in the laboratory. Experimental facilities 
designed to perform elastic electron scattering off exotic nuclei are at present under construction 
at RIKEN in Japan and GSI in Germany D24j-l2VD - 



2. Muonic atoms 



Muonic atoms were actually the first to provide accurate measurements of nuclear size [28], at a 
time when the nuclear radius was still believed to be 1.4 A 1 / 3 . The muon is a heavy unstable lepton 
with negative charge — e, and mass which is about 207 times heavier than that of the electron. A 
muon can form a bound system with a nucleus, just like the electron. However, because of its much 
heavier mass, the muon comes much closer to the nucleus, acting as a probe — in the Is state, for 
instance, the muon spends a significant fraction of time inside the nucleus. The ionization potential 
and the spectral lines corresponding to transitions between bound states are also much higher than 
the ones encountered in electronic atoms. For nuclei with large Z, the ionization potential is on 
the order of several MeV, and the transitions between muonic bound states range between 7 rays 
and X rays. 

When the muon is stopped in matter, it first scatters from atom to atom just as a free electron, 
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until it is captured into a high atomic orbit (a Rydberg state) of a particular atom. From this excited 
capture state the muon will cascade down through the electron cloud, ejecting Auger electrons and 
emitting radiation, until it arrives to a region inside the innermost electron orbit. From there, with 
no further electrons in sight, the muon decays mostly radiatively. The energies of these transitions 
provide information about the size and shape of the nucleus as well as tests of our understanding 
of QED. The capture and cascade processes up to the muonic Is state take place in 10~ 12 — 10~ 9 s, 
while the lifetime of the muon is of 2.2x 10~ 6 s. The muon spends thus most of its lifetime in the 
Is state, probing the nucleus. 

A comparison of the theoretical and experimental values for the energies and relative rates 
for various radiative transitions provide information about the nuclear size and the RMS radii. 
The muonic wavefunctions can be usually obtained with good accuracy, since the muonic atom 
is basically a hydrogen-like system — there is always only one muon bound to a nucleus, and the 
presence of the atomic electrons far away from the nucleus can be safely neglected. However, 
unlike in electronic atoms, the muon orbit is comparable to the nuclear size, and the transition 
energies of the same order of magnitude as the nuclear excitation energies, such that the nuclear 
polarization effects are considerably large. Also, even though the nucleus is still considerably 
heavier than the muon, the effect of nuclear motion (nuclear recoil, described in subsection III Al) 
should be taken into account [4]. It is beyond the purpose of this review to present the theory of 
muonic atoms and all the contributions to the transitions energies. The interested reader is referred 
to the review of Borie and Rinker ||29(]. 

The procedure for deducing the nuclear charge radii from the measured transition energies 
assuming a Fermi charge distribution has been described in detail in Ref. 13011 . which also gives 
a compilation of RMS radii, isotope and isomer shifts and magnetic hyperfine constants from 
muonic atoms. A more recent compilation with more precise data (also using a combined analysis 
of all electromagnetic methods to determine nuclear charge parameters) can be found in Ref. Ill2n . 
It should be noted that the main accuracy for the RMS radii deduced from transition energies 
in heavy muonic atoms is not given by the experimental errors, but by the uncertainty of the 
calculated nuclear polarization contribution. 

Just as in the case of electron scattering, RMS radii from muonic atom experiments have been 
determined for all stable isotopes already many years ago. Radioactive muonic atoms nowadays 
open the possibility to address exotic nuclei far from the valley of stability, whose properties are 
yet to be investigated. New, intense muon sources, together with techniques related to cold films of 
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radioactive atoms as host material, merging of ion and muon beams, or trapping of exotic isotopes 
and muons in Penning traps, will make a number of interesting muonic atoms experiments possible 
(see Ref. [31] on muon physics possibilities at muon-neutrino factories and references therein). 



3. Optical isotope shifts 

Optical isotope shifts are mostly obtained from laser spectroscopy on atoms or ions. Laser 
spectroscopy involves transitions with optical frequencies, which means that (with the exception 
of very light nuclei, discussed separately) this time outer atomic electrons get to probe the nuclear 
structure. The high resolution obtained by using laser-atomic beam spectroscopy and the high sen- 
sitivity of optical detection techniques allow precision measurements of the optical isotope shifts 
and hyperfine structure of stable as well as of rare and radioactive isotopes. Laser spectroscopy 
can also provide an alternative isotope separation method [32]. The advantage of this experimental 
method is its precision and the fact that it allows the determination of isotope shifts also for unsta- 
ble nuclei. The drawback is that the data is more difficult to interpret due to the presence of many 
electrons. In order to extract the most interesting nuclear parameter, the RMS radius, one needs to 
understand and discern the mass and field shift contributions out of the total isotope shift. This be- 
comes very difficult for many-electron atoms with large mass shift, for which the theory does not 
offer reliable SMS predictions. Usually this problem is solved by performing a combined analysis 
of isotope shift results from different types of experiments, and extracting scaling factors by hand 
by means of King plots (two-dimensional diagrams that represent isotope shifts obtained from two 



different methods on the x and y axis [13311 for purposes of calibrating the mass and field shifts in 
optical spectra). Unfortunately, the combined analysis is not possible for radioactive isotopes, for 
which no isotope shifts were determined from electron scattering or muon atom experiments. 

The typical optical isotope shift measurement involves a tunable laser which excites the op- 
tical transition of interest in atoms or ions. The transition energy is the object of a very precise 
determination by monitoring the excitation as a function of the laser wavelength. Usually either 
fluorescence or resonant ionization and subsequent ion detection (so-called Resonance Ionization 
Mass Spectrometry, RIMS) in the case of laser spectroscopy of atoms are used to detect the in- 
duced excitation. The typical problem which affects the accuracy is the Doppler effect. The atoms 
(usually in a beam) will have a velocity spread which shifts the resonance energy for the laser 
photons and thus broadens the fluorescence line. Doppler-free measurements are achieved by two- 
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photon laser spectroscopy, using one collinear and one anticollinear laser beams [320. Indeed, if 
the two absorbed photons with equal frequency u\ = 0J2 = oj, travel in opposite directions (their 
wave vectors are kx = —k 2 ), all particles, independent of their velocities, absorb at the same sum 
frequency u\ + u 2 = 2u. 

One unusual measurement of optical isotope shifts that does not involve laser spectroscopy 
has been already mentioned in Section III Al The experiment was performed on B-like and Be- 
like argon ions and confirmed the necessity of considering the relativistic nuclear recoil for the 
mass shift theory. The isotope shifts of 36 Ar versus 40 Ar were determined by measuring the 
ls 2 2s 2 2p 2 Pi/ 2 - 2 P-s/2 transition in Ar 13+ and the ls 2 2s2p 3 P 1 - 3 P 2 transition in Ar 14+ for the 
two isotopes with sub-ppm accuracy [8]. These optical frequency transitions both have magnetic 
dipole (Ml) multipolarity, such that the excited states, populated by collisions with the electron 
beam, are long-lived (or metastable) states. The experiment was performed at the Electron Beam 
Ion Trap (EBIT) in Heidelberg, Germany. The recorded spectrum is shown in Figure [2j A dedi- 
cated theoretical calculation showed that the NMS and SMS and their relativistic corrections are on 
the same order of magnitude. Furthermore, the theoretical isotope shifts including the relativistic 
corrections showed an excellent agreement with the experimental data, confirming the necessity to 
consider the nuclear recoil relativistically. Recently, the wavelength of the ls 2 2s 2 2p 2 P\/2 — 2 Pz/2 
transition in Ar 13+ has been confirmed by the same group in Heidelberg, this time using laser 
spectroscopy. For a comparison, the determined fluorescence spectrum for the it^, 2 Zeeman com- 
ponents is shown in Figure |3j 

A special case of optical isotope shifts are the very light nuclei (Z < 10), for which the field 
shift represents only a very small fraction of the total isotope shift (for Li, for instance, this fraction 
is 10~ 4 ). This means that RMS radii can be extracted from optical isotope shift measurements 
only if the accuracy of the experiment and the accuracy of the calculated mass shifts are on the 
order of 1CT 5 or better. Until recently, the SMS could not be calculated with sufficient accuracy 
for atomic systems with more than two electrons. For this reason, RMS radii determination for 
light atoms was restricted to electron scattering or muonic atoms experiments, which only work 
for stable isotopes. Laser spectroscopy experiments could provide RMS radii only for hydrogen 
and helium isotopes, for which the SMS was calculated with the necessary accuracy. The proton 
radius, for instance, has been determined from the Is — 2s transition frequency in hydrogen with 
10~ 14 accuracy D4I1 in the group of T. Hansch in Munich, Germany. 

Sufficiently accurate calculations of mass shifts for systems with up to three electrons became 
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reality recently Il35l - l39h . These calculations have enabled also the study of halo nuclei, a special 
variety of light nuclei that has received great attention since their discovery in 1985. In the n Li 
isotope it was found that the matter radius is much larger than that of its neighbours, suggesting 
a large deformation and/or a long tail in the matter distribution Il40ll . Later, such isotopes became 
known as halo nuclei, which have diffuse outer neutron distributions at the limit of stability. In 
particular, n Li is a three-body system with a 9 Li core and two halo neutrons, as depicted in Fig- 
ure HI Nuclear forces are not strong enough to bind one neutron to the 9 Li nucleus, nor are they 
able to bind two neutrons together in a dineutron. Yet the 9 Li and the two neutrons together can 
form a bound nucleus, n Li, with lifetime r = 8.4 ms. 

A feature of the halo nuclei is that the matter radius (given by the two halo neutrons) and 
the charge radius (of the three protons concentrated in the 9 Li core) do not coincide. Matter 
radii can be deduced from nuclear reactions. High-precision isotope shift measurements deliver 
the RMS charge radius. These two quantities are most important for understanding details of the 
nuclear structure and halo-core interactions. The RMS charge radius of 8 Li, 9 Li and n Li have been 
obtained by high-resolution laser spectroscopy on an atomic beam using two-photon Doppler-free 



excitation and resonance-ionization detection |41 , 



42]. Other laser spectroscopy 



the determination of the charge radii of other halo nuclei such as He 114311 - Be [44] or Ne [45] 



experiments for 



in ion beams or traps have been also performed. Ab initio descriptions of the structure of light 
nuclei based on realistic two- and three-nucleon interactions between individual nucleons have 
been developed during the past years, with work still in progress. 

4. Isotope shifts from x-ray transitions 

Isotope shifts from x-ray transitions are provided by x-ray spectroscopy experiments with atoms 
or highly charged ions. The K a x-ray transitions are the most obvious candidates. The very first 
such experiment in 1965 determined the isotope shift in K a fluorescence of 233 U and 238 U neutral 



atoms which were irradiated by broadband x-rays 114 611 . While from the theoretical point of view 
isotope shifts of atomic Ka x-ray transitions are much easier to handle, the experimental precision 
of x-ray spectroscopy is considerably lower than the one in the optical regime, where lasers are 
available. With the operation of the x-ray free electron laser LCLS at SLAC in Stanford, USA, 
providing photons with up to 8 keV energy, high-precision experiments in this energy regime will 
hopefully become possible in the near future. 
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Compared to experiments with neutral atoms, precision measurements of x-ray transitions in 
highly charged ions are even more advantageous due to the substantially simplified theory to ex- 
tract the nuclear charge parameters from the isotope shifts. X-ray spectroscopy in highly charged 
ions of different uranium isotopes has been performed at the SuperEBIT at the Lawrence Liv- 
ermore National Laboratory in the USA Ii47ll . The ion beam of the EBIT was used to ionize, 
excite and radially trap the ions, and the characteristic x-rays of the 2s 1/2 — 2p 3 / 2 transition 
of nearly 4.5 keV for Li-like, Be-like, B-like, C-like and O-like 233 U and 238 U ions were de- 
tected with a curved-crystal spectrometer with an energy resolution of 1 . 1 eV. This resolution of 
100 ppm in highly charged ions was already good enough to improve the accuracy of the RMS 
radii for uranium by about one order or magnitude. The obtained difference in the RMS radius 
#(r 2 ) 233 ' 238 = —0.457 ± 0.043 fm 2 lies in between earlier measurements based on different tech- 
niques. 



5. Dielectronic recombination experiments 

Recently, another type of atomic physics experiments involving electron recombination in 
highly charged ions have been used to obtain high-precision isotope shifts. This method involves 
a resonant channel of photorecombination, namely, dielectronic recombination. When a free elec- 
tron recombines into an ion, it can release its excess energy and angular momentum through a 
photon. This is the direct process of radiative recombination (RR), possible at all incoming elec- 
tron energies. However, there also exists a resonant channel: free electrons with matching kinetic 
energy can recombine resonantly and excite a bound electron, as shown schematically in Figure [51 
The formed excited electronic state can decay either by an Auger process (which is the reverse of 
the dielectronic capture), or by emitting one or more photons. The dielectronic capture followed 
by the emission of photon(s) from the bound electron(s) is known as dielectronic recombination 
(DR). DR is of great importance as the predominant electron recombination channel in solar and 
other dense astrophysical plasmas. Its study in highly charged ions provides important information 
about transition energies, isotope shifts, hyperfine splitting, in one word, electronic structure. By 
studying the electron recombination rate as a function of the free electron energy one obtains the 
spectra of DR resonances, whose positions give the atomic transition energies. 

DR experiments involve highly charged ions, either in a storage ring like the Experimental Stor- 
age Ring (ESR) at the GSI in Darmstadt, Germany, or in an EBIT. Depending on the experiment, 
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the photons following the recombination or the recombined ions are detected. In any case, both 
RR and DR channels are observed, with the resonance structure of DR on top of the fiat RR depen- 
dence on the free electron energy. An example is presented in Figure [61 which shows the photon 
energy distribution following recombination of mercury ions ranging from He-like to Be-like in 
an EBIT DR experiment. The electron energy was scanned across the region of the KLL reso- 
nances. From their positions, the atomic transition energies can be determined. In a storage ring 
recombination setup, the free electrons for DR are provided by a co-propagating electron target or 
by the electron cooler used to cool the ion beam. Particularly at very low ion-electron collision 
energies, storage ring experiments have enabled precise determination of Lamb shifts in several 
highly charged ions [48]. In order to determine very precise isotope shifts, DR of an electron with 
very small kinetic energy, preferably less than 10 eV, is required. For some nuclear charges Z, this 
is possible by allowing the capture into a Rydberg state to match a core n = 2 intrashell transition. 
Recombination into Rydberg states also has the advantage that the loosely bound Rydberg electron 
itself is unaffected by the small variations of nuclear potential due to different charge distributions 
and its contribution to the isotope shift can be safely neglected. 

Recently, isotope shift measurements in a DR experiment and subsequent RMS radii determi- 
nation have been achieved, opening a new direction for isotope shift experiments. Isotope shifts 
have been measured using the transitions e~ + Nd 57+ (ls 2 2si /2 ) -> Nd 56+ (ls 2 2p nlj)** of two 
Nd isotopes 142 Nd and 150 Nd |49|]. Here ** denotes the doubly excited state created in the dielec- 
tronic capture. The resonant electron capture occurs in a Rydberg state with principal quantum 
numbers n > 18 for the 2s — > 2pi/ 2 excitation and n > 8 for the 2s — > 2p 3 / 2 excitation. Figure [7] 
shows the dielectronic recombination rates for the two isotopes as a function of the recombina- 
tion energy, with visible shifts of the resonance positions. The interpretation of these shifts in 
terms of RMS radii is in this simple three-electron system very reliable and could be done within 
state-of-the-art structure calculations that take into account relativistic and QED contributions as- 
suming a two-parameter Fermi charge distribution for the nucleus. Out of the measured isotope 
shifts for the 2s - 2p 1/2 (5£ 142 - 150 =40.2(3)(6) meV) and 2s - 2p 3/2 (5£ 142 - 150 =42.3(12)(20) 
meV) transitions, with statistical and systematical error bars given in the first and second bracket, 
respectively, the RMS radii difference was determined, ^^ r 2 ^ 142 - 150 = — 1.36(1)(3) fm 2 . This 
value is slightly larger than the average experimental value of —1.291(6) fm 2 from the combined 
analysis in Ref. [19], but still in good agreement with the upper end of the individual data that 
build the basis of the evaluation. 
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III. HYPERFINE SPLITTING 



Hyperfine splitting (HFS) of atomic levels has been first observed by Michelson in 1881, but it 
was first in 1924 that Pauli suggested the observed effects might be due to magnetic interactions 
between electron and nucleus, possibly possessing a total angular momentum or "nuclear spin". 
Of course, at the time the concept of "spin" was not really established in quantum mechanics, but 
analogously to the fine structure that was attributed to the electron orbital angular momentum, the 
hyperfine features could be explained by assuming the existence of a nuclear magnetic moment. 

The HFS effects are due to the interaction of the nuclear multipole moments with the electro- 
magnetic field created by the electrons at the nucleus. The most important of these moments are 
the magnetic dipole moment, associated with the nuclear spin itself, and the electric quadrupole 
moment, which gives a measure of the deviation from a spherical charge distribution in the nu- 
cleus, i.e., of the nuclear deformation. For an atom or ion with nuclear spin / and electronic total 
angular momentum J, the value of the energy level splitting is given by (see, for instance, Ref. yj]) 



AP A B lK(K+l)-2I(I + l)J(J+l) 

AE =2 K + J J(2/-1)J(2J-1) ' (16) 

where K is given by 

K = F(F + 1)-I(I + 1)-J(J + 1), (17) 

and F is the quantum number of the total angular momentum F = I + J. The first term in Eq. (fT6l) 
is the magnetic dipole HFS, which is responsible for the splitting of one fine structure atomic 
energy level into a HFS multiplet labeled by F. The number of hyperfine structure components 
corresponding to a fine structure energy level is the smaller of the two numbers (2J + 1) and 
(21 + 1). The hyperfine separation, i.e., the energy difference between two neighbouring hyperfine 
levels is then given by AF, with the quantity A in Eq. (fT6l) given by 



A = ^9iWs J{J+ l ){2l + l) -^ , (18) 

where fi is the magnetic constant or vacuum permeability, gj the nuclear g factor, fi B and ijln the 
Bohr and nuclear magneton, respectively, / the orbital quantum number, n the principal quantum 
number and a M the reduced Bohr radius, a M = a (m/yu) obtained by using the reduced mass fx of 
the electron with respect with the nucleus. From the number of lines and the size of the splitting 
the nuclear spin and nuclear magnetic dipole moment can be extracted. 

The second term in Eq. (fT6l) represents the electric quadrupole correction. Because the de- 
pendence on the quantum number F is different from that of the magnetic dipole correction, the 
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hyperfine separation will no longer be simply proportional to F. The quadrupole coupling constant 
B is proportional to the electric quadrupole moment Q and the average gradient of the electric field 
produced by an electron at the nucleus, 

d 2 V 

The electric quadrupole contribution provides information about the spectroscopic quadrupole 
moments and nuclear deformation parameters. 

An illustrative example of the HFS scheme for 1 = 3/2 and J = 1 is presented in Figure [8] 
Since the hyperfine energy corrections are independent of the quantum number M F corresponding 
to the projection of the total angular momentum F on the quantization axis, the hyperfine levels 
show a (2F + l)-fold degeneracy which can be lifted in a magnetic field. Note that transitions 
between different HFS levels corresponding to the same fine structure level are forbidden in the 
dipole approximation. They occur as magnetic dipole (Ml) or electric quadrupole transitions 
(E2), and can be measured by radio-frequency methods. Due to the very good precision of such 
measurements, hyperfine transitions are used as time standards - the SI fundamental unit of time, 
the second, is defined in terms of the frequency of the transition between the levels with F = 4, 
M F = and F = 3, M F = of the ground state of 133 Cs. 



A. Isotope shifts in the hyperfine structure 

Isotope shifts were first observed in the HFS, and already in 1931 it was pointed out that the ob- 
served displacements were greater than expected from the mass correction of the Rydberg constant 
Jlp| . It was suggested that these shifts are due to deviations of the nuclear electric field from the 
Coulomb law. This represented the starting point for the whole theory of the field shift, outlined in 
Section HH For HFS, the discussion should be extended to the nuclear magnetization distribution 
and the effects of the extended nuclear charge distributions on magnetic interactions. 

In order to introduce the corrections to the HFS expression in Eq. (fT6l) . let us consider the 
simplest case of hydrogen-like ions, with one Is electron and a nuclear spin I. The relativistic 



magnetic HFS for the sublevel with largest F can be written as [|50D 



. _ — . a 4 Z 3 ujm e (I + J)m e c 2 r 4/ rWl s , 

AE = TT^ ^DPti) 1 ^" -fld-«) + *J • ao) 

where a is the fine structure constant, /ij the nuclear magnetic moment, A(ctZ) a relativistic 
correction factor and n iad QED radiative corrections. The product in front of the square bracket 
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is the AK/2 term from the non-relativistic expression in Eq. (fT6l) . evaluated for F = I + J. 
The remaining corrections 5 and e stand for the nuclear charge distribution and magnetization 
distribution corrections, respectively. 



In 1932, G. Breit and J. Rosenthal 115 111 studied the effect of the deviations of the nuclear field 
from the Coulomb law on the interaction of the electron with the nuclear spin. In this calculation, 
the nuclear magnetization was assumed to be concentrated at the center of the nucleus. The ra- 
tio between the HFS interaction assuming a realistic nuclear charge distribution and the one with 
a point-like nucleus is known as the Breit-Rosenthal correction 5, sometimes simply called the 
nuclear charge distribution correction for HFS. The difference arises from the change of the elec- 
tronic wave functions when assuming a nuclear charge distribution instead of a point-like charge, 
just as in the case of the field shift. 



In 1950, A. Bohr and V. Weisskopf 115211 studied the influence of nuclear structure on the HFS 
of heavy elements. The issue started from the so-called "hyperfine anomaly", an unexplained 
isotopic dependence in the ratio between HFS and nuclear g factors. For the first time they took 
into account the possibility that the nuclear dipole moment is not point-like, but should be repre- 
sented by some distribution of magnetism over the nuclear volume. Indeed, since the total nuclear 
spin results from the individual spins of neutrons and protons and their angular momenta, the 
distribution of magnetic dipole density over the nuclear volume may vary greatly from nucleus 
to nucleus depending on the relative spin and orbital contributions to the total nuclear moment. 
The ratio between the HFS interaction energy assuming the extended nuclear magnetization and 
the one for point-like nuclear dipole moments is known as the Bohr- Weisskopf correction e. For 
heavy H-like systems the Breit-Rosenthal and Bohr- Weisskopf corrections are much larger than 
the QED corrections. For the case of H-like 165 Ho, for instance, vacuum polarization and self- 
energy contributions (included in K ra d) account for 0.5%, the nuclear magnetization distribution 
(Bohr- Weisskopf) for 2% and the Breit-Rosenthal effect for 8% of the total energy splitting. 

Accurate HFS measurements in H-like systems provide a way to determine nuclear magnetiza- 
tion distribution radii. By using previously determined nuclear magnetic moments, and applying 
appropriate corrections for the nuclear charge distribution and radiative effects, the experimen- 
tally determined HFS can be interpreted in terms of nuclear magnetization radii. The aim is to 
distinguish between models of the magnetic structure of the nucleus. This is another opportu- 
nity to learn more about the nucleus via high-precision atomic physics experiments assisted by an 
adequate theoretical interpretation. Since often several effects contribute to the atomic transition 
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energies, increased experimental and theoretical precision refines the search for other features at 
the borderline between atomic and nuclear physics. For instance, the uncertainty in the neutron 
distribution inside the nucleus affects the interpretation of parity non-conservation experiments, 
and knowledge about the nuclear magnetization distribution is vital for the interpretation of exper- 
iments looking for nuclear anapole moments or a proton electric dipole moment. 



B. Highlights of hyperfine splitting experiments 

The year 2007 marked the 40 th anniversary of the redefinition of the International System of 
Units (SI) second based on a hyperfine transition in cesium. Until 1967, the second had been 
always defined based on astronomical time scales. The hyperfine transition F = 4, Mp = 
and F = 3, M F = of the ground state of 133 Cs has a defined frequency of 9192.631770 GHz, 
commensurate with the accuracy in measuring this transition frequency achieved in microwave 
experiments. 

The original idea of microwave atomic beam magnetic resonance experiments (or, as they were 



called at the time, radiofrequency experiments) is the following p3LI54[) : a beam of atoms propa- 
gates through two inhomogeneous magnetic fields oriented such that the deflections they cause are 
equal but opposite. The net deflection is then zero if and only if the magnetic moment of the atom 
in the direction of each field is the same. If the orientation of the atomic magnetic moment changes 
in the space between the two magnets, the atoms will not reach the detector placed at the zero de- 
flection position, and the latter will register an intensity drop. Between these inhomogeneous 
fields a region of homogeneous variable magnetic field is placed and perpendicular to this, a weak 
radiofrequency oscillating field is present. If the oscillating frequency corresponds to the energy 
difference between two hyperfine levels of interest, detectable transitions changing the orientation 
of the atomic magnetic moment will occur. Thus, the beam intensity dependence as a function of 
the applied frequency provides the atomic radiofrequency spectrum. A further improvement by 
Ramsey [55] utilizes two separated oscillatory fields basically building up an interferometer for 
the atomic beams. Nowadays, the oscillating field is applied in a microwave interrogation cavity 
which allows narrowing of the resonance line, reduces the sensitivity to microwave power fluctua- 
tions and magnetic fields by orders or magnitude, and it eliminates the Doppler effect. The Ramsey 
method yields the until recently unrivalled accuracy of the cesium beam primary frequency stan- 
dards, which achieved with NIST-7 a fractional uncertainty of 5 x 10~ 15 [56]. 
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Microwave experiments were the first ones to drive transitions between hyperfine atomic levels. 
Apart of cesium and other alkali atoms suitable because of their single valence electron, in particu- 
lar the 1 .4 GHz splitting frequency in hydrogen has been intensively studied, eventually providing 
the hydrogen maser with a relative frequency accuracy of 7 x 1CT 13 [57]. This illustrates the po- 
tential for extending ground state HFS investigation to H-like ions, due to the simplicity of the 
theoretical modeling for one-electron systems. However, the microwave measurement of the Is 
splitting in 3 He + 115 811 has not been followed by further steps towards higher- Z ions, mainly due 
to difficulties in producing and storing highly charged ions. Obviously, due to the energy splitting 
dependence on Z 3 , the HFS transitions in highly charged ions are also no longer accessible to mi- 
crowave spectroscopy. At around Z ~ 65, the transition between the two F levels of the Is ground 
state already scales up into the infrared, the visible and UV region of the electromagnetic spec- 
trum, depending on the particular value of the magnetic moment of the nucleus. High-precision 
laser spectroscopy then becomes the appropriate tool to address such transitions. 

Similar to the case of isotope shift measurements, HFS have been determined also in muonic 
atoms or in x-ray atomic spectra. Ground state HFS has been observed in x-ray spectra recorded 
after nuclear electron capture from the K-shell of cesium atoms [59]. The precision of this mea- 
surement was not sufficient to draw any conclusions about the nuclear magnetization. For a long 
time, the main tool of studying the nuclear magnetization distribution has been by measuring 7 
rays emitted in muonic atoms I160ll . The HFS of the ground state is in muonic atoms on the order 
of a few keV, which unfortunately have to be extracted from 7-ray transitions corresponding to the 
nuclear or muonic transitions, with energies of several MeV. The ground state HFS usually cannot 
be resolved directly, but it has to be determined from the broadening of the line using calibration 
lines of isotopes with no nuclear spin Q60Q . The HFS of the Is level can be determined with sev- 
eral percent accuracy, and since the effect of the nuclear magnetization distribution reduces the 
separation by 50 %, a value accurate to 10 % can be obtained for the Bohr-Weisskopf correction. 
The ground state HFS have been measured in muonic atoms to a precision of a few percent in 
209 Bi [|6l|l and 203 > 205 T1 [|62|l . Unfortunately, the uncertainties do not allow to distinguish between 
different nuclear models. 

For electronic atoms, the advantage of determining HFS in H-like ions is obvious. The HFS 
is large, the Is electron is close to the nucleus and feels the nuclear magnetization distribution 
strongly, and the theoretical interpretation of the measured spectra is reliable. Furthermore, nu- 
clear polarization effects are smaller than in muonic atoms. As soon as H-like heavy ions could 
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be produced, the two types of highly charged ions facilities, the storage ring and the EBIT, aimed 
at performing HFS experiments. Out of this competition, HFS for the ground state of 165 Ho 66+ 

fl, 185 < 187 Re 74+ Hi, 203,205 T1 80+ [65L 207p b 81+ and 209^82+ gj ^ determined The 

experiments at the ESR at GSI Darmstadt, which determined the HSF of 209 Bi 82+ (visible) and 

207p b 81+ ( infrared ) 

were performed by precision laser spectroscopy on the electron-cooled beam 



of H-like ions in a collinear antiparallel setup 11661 . 



67|]. Laser fluorescence was detected as the 



laser frequency was tuned in and out of resonance with the Doppler-shifted transition, with total 
accuracy limited by ion velocity uncertainties. Note that a certain a priori knowledge of the hy- 
perfine transition is required, since otherwise the scan in the laser frequency is prohibitively time 
consuming. 

At the EBIT there is no such impediment because laser spectroscopy is not needed. The hyper- 
fine transition of the Is ground state level is determined using "passive" emission spectroscopy, 
since the transitions are excited by electron collisions with the electron beam itself. To match the 
precision obtained by laser spectroscopy on highly charged ions at the GSI, a sensitive spectrome- 
ter with good resolving power was used. In Figure |9] the measured signal in the direct observation 
of the spontaneous emission of the hyperfine transitions F = 4 to F = 3 in the ground state of hy- 
drogenlike 165 Ho 66+ at the Livermore EBIT Q63D is shown. The HFS in 165 Ho 66+ , 2 °3,205 T1 80+ and 
185 ' 187 Re 74+ were determined at the Livermore EBIT with accuracies comparable or better than 
those of the laser-fluorescence measurements of 209 Bi 82+ and 207 Pb 81+ at the GSI. The nuclear 
magnetization radii deduced from EBIT experiments for the rhenium and thallium isotopes can be 



found in Refs. Q64 , 



650. 



Last but not least, DR can also be used to determine HFS of highly charged ions. The precision 
required by the HFS experiments is usually not matched by the typical accuracies of DR electronic 
transition energies. However, as already discussed in Section HI~Cl precision can be enhanced using 
resonances that lie very close to the ionization threshold. With relative velocities between ion and 
recombining electron very close to zero, the accuracy of the free electron energy is substantially 
increased. While there is no possibility to obtain HFS for H-like ions in DR, which always involves 
two electrons, experiments with other isoelectronic sequences offer the possibility to obtain precise 
transition energies and resolve HFS for different ion charges where QED, nuclear magnetization 
distribution and nuclear charge distribution contribute differently. 

Determination of the HFS of 4si/ 2 and 4pi/ 2 in 207 Pb 53+ [|68ll was possible by comparing DR 
resonances of 207 Pb 53+ and 208 Pb 53+ , since the latter has zero nuclear spin and does not present 
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HFS. The experiments were performed at the heavy-ion storage ring CRYRING in Stockholm, 
Sweden. The comparison between the observed resonances in the two lead isotopes was per- 
formed by taking into account all isotope dependent contributions using relativistic many-body 
perturbation theory and a multi-configuration Dirac-Fock method 116 811 . A careful evaluation ex- 
tracted the HFS with an accuracy on the order of 10% of the HFS constant. 

DR close to the recombination threshold can actually be used as a very accurate tool to obtain 
electronic transitions energies. Recently, DR Rydberg resonances below 0.07 eV were determined 
for 45 Sc 18+ with absolute accuracies below 0.0002 eV, a factor of 10 smaller than previous storage 
ring studies [69]. This remarkable precision (4.6 ppm) in the measured value of the 2si/ 2 — 2p 3 / 2 
transition energy enabled features of the HFS of the 2s state to be resolved. In conjunction with 
accurate atomic structure calculations, this experiment was sensitive to QED screening corrections. 

Finally, DR has been also used to determine the hyperfine-induced 2s2p 3 P — >■ 2s 2 1 So tran- 
sition rate in divalent beryllium- like Ti ions [70] at the heavy-ion storage ring TSR in Heidelberg. 
Hyperfine-induced transitions are forbidden electronic transitions rendered possible by the mag- 
netic interaction with the nuclear spin. For instance, it is well known that one-photon transitions 
between two states for which the electronic wave functions both have a total angular momentum 
equal to zero are forbidden by selection rules of any electric or magnetic multipole order. These 
transitions can occur, with very low probabilities, via two-photon emission. However, if the nu- 
cleus has a spin I different from zero, the hyperfine interaction mixes states with different angular 
momenta and thus one-photon transitions are no longer forbidden. Such hyperfine induced tran- 
sitions are interesting for obtaining ultra precise optical frequency standards and for cold-atom 
studies, or in astrophysics, where they can be used for diagnostics of low-density plasmas. 

The rate of the hyperfine-induced 2s2p 3 P — > 2s 2 1 S in 47 Ti 18+ was determined by comparing 
DR spectra of two Ti isotopes, 47 Ti 18+ and 48 Ti 18+ , with nuclear spins I = 5/2 and 1 = 0, 
respectively. The ions were produced mainly in the ground state 2s 2 1 S , with about 5% in the 
metastable 2s2p 3 P state. The electron-ion collision energy in the storage ring was tuned to a 
value where DR occurs only for ions in the excited metastable 2s2p 3 P state, and the rate of 
recombined ions Ti 17+ was monitored as a function of the storage time. The ions were stored in 
the ring for about 200 s. During this time the two-photon decay of the metastable state, with an 
estimated lifetime of approximately 4 days 117 ill , did not play any role. 

While for the 48 Ti ions the — > transition is forbidden, for the 47 Ti ions the lifetime of 
the metastable state is substantially reduced and the DR resonance essentially disappears. Out of 
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the recombination rates as a function of the ion storage time the hyperfine-induced decay rate of 
1.8 s was determined - about 40% lower than the theoretical value of 2.8 s. This determination 
is almost one order of magnitude more precise than the previous experimental value for the same 
2s2p 3 P — > 2s 2 1 S transition in N 3+ that was obtained from astrophysical observations and 
modeling [17211 . Hyperfine quenching of forbidden transitions has also been proposed as a method 
to probe the nuclear state. In collisions of 238 U, which has ground state nuclear spin equal to zero, 
a fast decay of the 2s2p 3 P state could provide evidence of the nucleus being in an excited state 
1731]. HFS experiments on excited nuclear states would be particularly interesting for identification 
of certain long-lived nuclear states, as pointed out later in Section |V] 



IV. EXOTIC NUCLEAR PROPERTIES: NUCLEAR EXCITED STATES AND THE WEAK IN- 
TERACTION 

So far, we have addressed changes in atomic transitions that are due to the nuclear size, mass 
and spin. However, these are not the only nuclear properties that can play a role in the life of an 
atom. More exotic features such as the nuclear excited states that determine nuclear polarization 
corrections for the atomic transitions, or of which excitation and decay mechanisms can directly 
involve atomic electrons, may have an important contribution to atomic transitions. Going to the 
weak interaction between nuclear constituents, the coupling of nucleons to the atomic electrons 
gives rise to atomic parity violation effects. 



A. Atomic parity violation 

In the absence of any external field to define a specific orientation in space, atoms do not have 
any preference whether to interact with left or right handed circularly polarized photons. Atomic 
states have a definite parity, and photon absorption or emission is equally probable for both left 
handed or right handed circularly polarized photons. 

This holds, however, only if the weak interaction between the atomic electrons and the nuclear 
constituents is neglected. The weak interaction can mix atomic states of different parity, producing 
observable effects: left and right handed photons are absorbed and emitted slightly differently, i.e., 
the atoms show optical activity, so-called circular dichroism. Parity is no longer conserved — one 
refers to parity violation or parity non-conservation in atoms. Such results are in conflict with 
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QED, and support the theory of unification of the electromagnetic and weak interactions. 

In quantum mechanics the parity operator P transforms each point of space f into its opposite 
— r. The operator P thus performs a space inversion and transforms a quantum state into its mirror 
image. A parity quantum number is associated to P, namely, p — ±1, given for atomic states by 
p = (—1)', where I is the orbital quantum number. The parity quantum number is a multiplicative 
quantum number, for instance, for a system of N independent electrons, it will be the product 
of the N individual quantum numbers pj = (— . For electromagnetic interactions, parity is 
conserved, meaning that in any transition the products of the parity quantum numbers of the initial 
and final sub-systems are identical. 

For a long time it was believed that parity conservation is a universal law, since there were 
no experiments to prove the contrary. In the early 1950s, a puzzling particle physics experiment 
was the first to question the universality of parity conservation. The K meson was found to de- 
cay sometimes in two pions and sometimes in three pions, and due to the p = — 1 parity of the 
individual pions, this implied that the two observed final states had different parity! Lee and Yang 
stepped into the unknown and suggested that parity is not conserved in the decay - since they 
noted that at the time there was no experimental evidence for parity conservation in transitions 
induced by the weak interaction. One year later, in 1957, Wu confirmed experimentally that parity 
is violated in the nuclear decay 117 511 . 

Shortly afterwards, the idea to search for parity violation in atomic transitions was considered 
H76J] . However, a first estimate of the left-right asymmetry in atomic transitions was so exceedingly 
small (1CT 15 ) that its observation seemed completely hopeless. The rotation of the polarization 
plane of visible light propagating through optically inactive matter was estimated to be only 1CT 13 
rad/m. Luckily, this estimate proved to be by far too pessimistic, and atomic parity violation was 
experimentally observed some decades later. 

Parity violation in the weak interaction was first observed in nuclear (3 decay, where always a 
change in the electric charge and consequently in the identity of the decaying particles occurs. It 
was therefore difficult to conceive that the weak interaction and correspondingly its parity viola- 
tion feature could have any influence in neutral stable atoms, where the identity of the interacting 
particles is obviously preserved. Zel'dovich Ir76ll . when giving the first discussion on the pos- 
sible effects in atomic physics of weak electron-nucleon interaction, referred to so-called neutral 
currents. 

It took many years for the theoretical understanding of the weak interaction to become math- 
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ematically satisfactory — by the late 1960s, Glashow [TZO, Weinberg 117 811 and Salam and Ward 



[79] had suggested independently that the weak and the electromagnetic interactions can be un- 
derstood as different manifestations of a unified underlying electroweak interaction. The family 
of gauge bosons, the particles that mediate the weak interaction, was known to contain the W + 
and W~ bosons, both carrying a unit of electric charge. The W ± gauge bosons mediate the weak 
interaction in (3 decay. With the theoretical unification of the weak and electromagnetic interac- 
tions, an important prediction was the existence of a new gauge boson, called Z°, mediating weak 
neutral current interactions. The Z° boson carries no electric charge, so one way it couples to 
the fundamental particles such as electrons and quarks is like a heavy photon, as depicted in the 
Feynman diagram in Figure \\0\ However, an important aspect of the Z° boson is that it has chiral 
behaviour — it changes its coupling to constituents when reflected in a mirror — similar to its kins 
the W ± bosons. The coupling of Z° to the electron is proportional to the electron helicity h e , a 
quantity odd under space reflection, defined as the scalar product between the electron spin and 
velocity, h e — a e ■ v e . The short-range electron-nucleus Z° exchange potential thus leads to the 
parity violation effects in atoms. 

The Coulomb interaction and Z° boson exchange electron-nuclear potentials have certain sim- 
ilarities. In the vicinity of the nucleus (assumed to be pointlike), the Coulomb potential can be 
expressed as 

Ze 2 

V c (r e ) = — , (21) 



and the parity-violating potential can be written in the non-relativistic limit as [18011 



V P v{r e ) = -zQw9z° (c e • v e )/c + h.c. , (22) 

where g z o is a coupling constant whose size, as a consequence of the electroweak unification, is of 
the order of e. The strengths of the Coulomb and electroweak interactions between the electron and 
the nucleus are given by the Ze 2 and Qw9 z o factors, respectively. The quantity Qw is therefore 
known as the weak charge of the nucleus — in analogy to the nuclear electric charge. Due to the 
very small Compton wavelength h/(M z oc) of the Z° boson, it is legitimate to take the limit of 
an infinite Z° mass, which then allows replacing the Yukawa potential in Eq. (l22l) by a Dirac 
distribution, 

lim (M^) 2e M-M ZO cr e /h) = 

M z0 -+oo \ n J r e 
Usually the Z° boson exchange potential is written with the help of the Fermi constant, Gp = 
4 x lCT 14 Ryd x a^, where Ryd is the atomic Rydberg constant and a the Bohr radius. The 
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parity-violating potential can then be written as 

QwGf r r-3/ 



V pv (r e ) = ^-^[5\r e )(a e -v e )/c + h. C .]. (24) 



The Coulomb potential in Eq. (|2~TT) is even under space reflection, while the Z° boson exchange 
potential in Eq. (|24l) is odd, giving rise to left-right asymmetry in atomic transitions. 

The basic principle of most parity violation experiments is to compare the transition rate be- 
tween two states A and B with the one between their mirrored states A and B. The left-right 
asymmetry is defined then as the difference between these rates divided by their sum. With an 
even amplitude coming from the electromagnetic interaction and the weak amplitude that has an 
odd contribution, we have the transition rates Pl/r = \A em ± v4^ d | 2 , which leads to the left-right 
asymmetry 

A LR = ^— ^ = 2Re{A°4 d /A em ) . (25) 

While at first Alr was estimated at 1CT 15 , which seemed to make observation impossible, Marie- 
Anne and Claude Bouchiat 118 IN predicted that the electroweak effects grow slightly faster than the 
cube of the atomic number Z. Thus, heavy atoms are most appropriate to look for atomic parity 
violation, e.g., for cesium (Z=55) the parity violation effects are increased by a factor of 10 6 . The 
cesium atom still represents the best compromise between high Z and simple atomic structure 
allowing for precise atomic structure calculations. Another way to enhance the effect is the study 
of forbidden atomic transitions such as the 6sx/ 2 — 7si/ 2 transition in Cs atoms. This transition 
is a magnetic dipole, strongly suppressed because of the different principal quantum numbers of 
the two states. Parity mixing due to the weak interaction induces a small electric dipole (El) 
transition amplitude, which has been observed experimentally. The left-right asymmetry in this 
case can reach 10~ 4 . 

Unfortunately, due to the very small oscillator strength of this 6sj / 2 — 7si/ 2 transition, there were 
a number of practical problems to see any signal. The Bouchiats [82] suggested a clever technique 
for enhancing the signal relying on the interference of a Stark-induced mixing of opposite parity 
states in an atom and the mixing caused by weak neutral currents. The Stark-induced amplitude 
works as a lever arm to magnify the neutral current amplitude. Following this proposal, a number 
of groups endeavored to search for weak neutral currents in atomic parity violation experiments, 
using either the Stark-mixing idea or by studying the rotation of plane-polarized light as it passes 
through a gas of atoms. The reader is referred to the comprehensive review in Ref. 
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The most precise atomic parity-non-conservation experiment (experimental accuracy of less 



than 0.4%) Il83h examined the mixing of S and P states in atomic cesium. Specifically following 
the original proposal in Ref. [82], it compared the mixing due to the parity-violating neutral weak 
current interaction to the S-P mixing caused by an applied electric field. The measurements were 
interpreted in terms of the nuclear weak charge Qw, quantifying the strength of the electroweak 
coupling between atomic electrons and quarks of the nucleus. The main goal is to compare this 
experimental result with the lowest order predictions of the Standard Model for the weak charge, 
where Qw is the sum of all the weak charges of the atomic nucleus constituents, the u and d 
quarks, 

Q w = (2Z + N)Q w (u) + {Z + 2N)Q w {d). (26) 
The weak charge happens to lie close to the neutron number N, 

Qw = —N — Z(A sin 2 8 W — 1) — —N , (27) 

where the weak mixing angle Qw is a free parameter of the theory, determined experimentally to 
be sin 2 6 W = 0.23. 

In order to obtain the experimental value of Qw, atomic structure calculations for the amount 
of Stark mixing and the relevant parity-violating electronic matrix elements are needed. The un- 
certainties in the determination of Qw are actually dominated by the uncertainties in these two cal- 



culated quantities. In Ref. 118411 the evaluation of the experimental data in Ref. [18311 was improved 
by also determining experimentally the Stark mixing and incorporating further experimental re- 
sults into the evaluation of the uncertainty in the calculation of the relevant parity-violating matrix 
elements. The most accurate to-date determination of this coupling strength by combining these 
two measurements using high-precision calculations for the cesium atom yield a weak charge of 



Qw = — 73.16(29) ex p(20)theor B85D, which is in excellent agreement with the predictions of the 
Standard Model. 

Atomic parity-non-conservation measurements are uniquely sensitive to a variety of new 
physics, such as the existence of additional Z bosons, because they probe a different set of model- 
independent electron-quark coupling constants than those measured by high-energy experiments, 
in a different energy regime Il86n . Given the huge effort made in testing the Standard Model in 
high-energy collider experiments, the results of the atomic parity violation measurements repre- 
sent a significant triumph for table-top physics. Improved precision of the atomic experiments, 
possibly to be achieved by measurements for several atoms along an isotopic chain, as a way to 
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circumvent the atomic theory uncertainties [87|], is expected to keep parity-violation experiments 
in atoms interesting and useful for the continuing quest for new physics beyond the Standard 
Model. 



B. Nuclear polarization effects 

With increased precision of atomic physics experiments and QED as a powerful tool to calcu- 
late radiative corrections up to high orders in a, one would be tempted to think that it is only a 
matter of time until excellent agreement between theory and experiment will be achieved. This is 
however not the case, since there are also other corrections that have to be accounted for. The nu- 
cleus itself is one of the main sources of undesirable interferences. At a certain level of accuracy, it 
is no longer enough to consider the nucleus as source of the classical external Coulomb field mod- 
eled by a classical charge density distribution. One should also take into account the interaction 
between the radiation field and the nucleus due to intrinsic nuclear dynamics. Due to the exchange 
of virtual photons, the nucleus can undergo virtual transitions to excited states, giving rise to en- 
ergy shifts of electronic states. These contributions are called nuclear polarization effects. The 
nuclear polarization contribution is always present in atomic transition energies and represents a 
natural limitation for any high-precision tests of QED. As already mentioned in Section HTCl nu- 
clear polarization is often the largest uncertainty factor when it comes to identifying the various 
contributions to the measured isotope shifts. Nonetheless, one can imagine reaching a level of 
experimental and theoretical accuracy which allows testing the influence of internal nuclear struc- 
ture as well as higher-order QED effects. Until then, the nuclear polarization corrections have 
to be calculated theoretically up to a certain level of accuracy and then used in the complicated 
procedure of extracting from the experimental data the atomic or nuclear quantity of interest. 

Some general remarks are due here. The nuclear polarization correction arises from virtual 
transitions as shown in the Feynman diagrams in Figure QT] Since in the nuclear polarization 
energy shift expression the energy difference E a — E n between the atomic and nuclear virtual 



transitions appears in the denominator (880, the closer the two energies are, the larger the nuclear 
polarization correction. This explains why the first systems for which nuclear polarization effects 
have been studied are the muonic atoms, for which the energy corrections can be as large as several 
keV. In atomic spectra, on the other hand, nuclear polarization corrections have been disregarded 
for a long time as being too small to be experimentally observed. The muon transition energies 
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are of the order of MeV, on the same scale as typical nuclear transitions. This does not hold for 
atomic transitions, where the transition energies are usually below or around 100 keV, thus much 
smaller than the nuclear excitation energies. Exceptionally, there are a number of low-lying rota- 
tional and vibrational bands, mostly in even-even nuclei, which can also reach values below 100 
keV. Some very strong nuclear transitions in the MeV range, called giant dipole resonances, also 
play an important role for nuclear polarization corrections, not because of their match in energy, 
but because of their large resonance strength. Unfortunately, since nuclei have a complex inner 
structure and their excited level schemes present more exceptions than rules, for any particular 
nucleus under consideration an adequate choice of the relevant excitations and nuclear models to 
describe them is necessary. One usually considers only those having strong excitation strengths, 
that are expected to contribute considerably. 

Collective nuclear models are the simplest alternative, describing collective nuclear rotational 
and vibrational states as well as giant dipole resonances. These phenomenological collective mod- 
els of the nucleus have as underlying physical picture the classical charged liquid drop, disregard- 
ing the interior structure (i.e. the existence of the individual nucleons) in favor of the picture of a 
homogeneous fluid-like nuclear matter. The model is applicable only if the size of the nucleon can 
be neglected with respect to the size of the nucleus as a whole, as it happens in the case of heavy 
nuclei. Furthermore, collective models are only appropriate to describe even nuclei. For low-lying 
single-particle excitations of a valence proton in odd nuclei, nuclear polarization calculations rely 
on the Weisskopf approximation j89|l and require a different approach to calculate the effective 
electron-nucleus interaction. 



A review of the theory of nuclear polarization in muonic atoms can be found in Ref. [|29D. A 
relativistic, field-theory approach to nuclear polarization calculations for electronic heavy atoms 
has been put forward by Mohr, Plunien and Soff [90]. With the introduction of an effective photon 
propagator with nuclear polarization insertions, the nuclear polarization corrections are treated on 
equal footage as radiative corrections to the electron energy, and appear in an effective self-energy. 
It is beyond the scope of this work to present the full theoretical treatment of nuclear polarization 
corrections to atomic transition energies. Instead, the reader is referred to the detailed description 
of the effective photon propagator technique in Ref. [88]. For very light atoms such as deuterium, 



92D. 



the theory of nuclear polarization has been developed in Refs. H91L 

Because of the unreliable nuclear part, any calculation of nuclear polarization shifts is inher- 
ently phenomenological and depends on the parameters of the nuclear model used. On average a 
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typical error of 25% is assumed, most of it coming from the uncertainty in the nuclear parameters, 
e.g., the nuclear transition energies and strengths. The choice of relevant nuclear states is another 
source of uncertainty, since among the neglected ones some important contributor could have been 
missed. 

The nuclear polarization corrections are usually small, on the order of 200 meV for the K shell 
electron in g2 8 U (compared to the 200 eV shift due to the extended nuclear size). An evaluation of 
the nuclear polarization energy shifts for the K shell of g| 8 U and ^SfPb 118 811 shows that for I^Pb, 
with the first excited nuclear state already as high as 2614 keV, the nuclear polarization shift is three 
orders of magnitude smaller than in g2 8 U [93], which has a low-lying rotational band with a number 
of excited nuclear levels below 600 keV, the first at 44 keV. Nuclear polarization corrections have 
been calculated mostly for systems for which the interpretation of experimental data required high 
accuracy in theoretical calculations and inclusion off all contributions in the atomic transition 
energy shifts. Among these are the nuclear polarization correction for the recently studied halo 
nucleus 11 Li ||94fl . or for various helium isotopes [|95|], the nuclear polarization correction to the 
bound-electron g-factor in heavy H-like ions or for QED tests in heavy, highly-charged ions 



[97, 



98]. 



Unfortunately, small as the nuclear polarization corrections are, they remain a theoretical chal- 
lenge and a source of uncertainty for atomic structure calculations, being the major accuracy lim- 
itation for the determination of important nuclear parameters such as the RMS radii from atomic 
transition energies. 



C. Nuclear excitation mechanisms involving atomic electrons 

In some cases, excited nuclear states can directly drive atomic transitions. It is well known that 
for low-lying nuclear excited states, a frequent and sometimes preferred decay channel is internal 
conversion (IC). Instead of emitting a 7 photon, the nucleus transfers its excess energy and angular 
momentum to a bound electron, which leaves the atom. IC and 7 decay are usually competing 
nuclear decay channels. Their ratio defines the internal conversion coefficient, a = A IC /A^, 
where A IC and A y are the IC and 7 decay rates, respectively. IC is predominant for small transition 
energies (< 100 keV) and is the primary decay channel for + — > + forbidden transitions whose 
energy does not allow decay via pair production. 

The inverse process of IC can occur only when ions and free electrons come across, i.e., in a 
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process of electron recombination. The two most important channels of electron recombination 
have been discussed when introducing DR: if the electron capture into the atomic bound state 
occurs with the emission of a photon, one speaks of radiative recombination. When the capture 
occurs in the presence of a bound electron, at the resonance energy that permits recombination 
with the simultaneous excitation of the bound electron in the ion, one refers to dielectronic re- 
combination. The inverse process of IC, nuclear excitation by electron capture (NEEC), is the 
nuclear physics analogue of DR, in which the place of the bound electron is taken by the nucleus. 
If the electronic and nuclear energy levels match, the recombination can take place with the si- 
multaneous excitation of the nucleus, as shown schematically in Figure [T21 The nuclear excited 
state decays then radiatively or can undergo internal conversion. In contrast to DR, NEEC can also 
occur in bare ions, as the presence of a bound electron is not required. 

IC and NEEC also have more exotic siblings, the bound internal conversion (BIC) and its 
inverse process nuclear excitation by electron transition (NEET). Bound IC is a resonant nuclear 
decay channel which may occur if the nuclear excitation energy is not enough to ionize the atomic 
electron, but can induce a transition between two bound electronic states. Nuclear excitation by 
electron transition is a simultaneous excitation of the nucleus during an atomic decay transition 
when the energies of the two transitions match. There are only very few atoms in which such a 
match of the atomic and nuclear transition energies exist, see for instance the list of candidates 
for NEET in Il99n . Although it is difficult to find systems that fulfill the resona nce c ondition, 
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BIC and NEET have already been observed experimentally in the year 2000 [100, 
NEEC together with BIC and NEET couple atomic and nuclear transitions and convert electronic 
orbital energy into nuclear energy and vice versa. They offer therefore the possibility to explore 
the spectral properties of heavy nuclei, such as nuclear transition energies and strengths, through 
atomic physics experiments. NEEC and NEET are also expected to allow the study of atomic 
vacancy effects on nuclear level population and lifetimes. 

That the atomic configuration can influence the nucleus has been long known. Maybe the most 
spectacular case arises when a bound (3 decay channel is opened in highly charged ions. In (3 
decay the nucleus of an atom emits an electron into the continuum. If the nuclear decay lacks the 
energy to emit the electron into the continuum, the (3 decay channel is closed. This can change 
dramatically if instead of an atom we consider a highly charged ion. The (3 decay electron is then 
emitted in a bound state, requiring less energy for the nuclear transition. The opening of the new 
bound (3 decay channel in highly charged ions considerably influences the lifetime of unstable 
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levels in nuclei H102L I103H . As an example, the lifetime of the ground state 75 7 Re decreases by 
more than nine orders of magnitude from 42 Gyr for the neutral atom to 32.9 yr for bare ions as a 
consequen ce of new bound f3 decay branches to the ground and excited states of the 7g 7 0s daughter 
[| 1 04L 1 1 OSh . The case of ^f/Re is particularly interesting in astrophysical context, since it severely 
affects the accuracy of the 187 Re- 187 0s cosmochronometer II 10511 . The motivation in investigating 



the behavior of nuclei in highly charged ions is thus related to nuclear astrophysics, studies of 
nuclear decay properties and tests of the relativistic description of atomic inner-shell processes. 
For highly charged ions, also NEEC is expected to play an important role for the nuclear lifetimes 
and the excited-state population in high-density astrophysical plasma environments. 

Historically, the NEEC recombination mechanism has been proposed for the first time for laser- 
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106 



107D, solid targets 11108 



? ollow ed by 



UOborfor 



assisted plasma environments in 1976 in Ref. [110611 . Whi 
a number of theoretical studies of the process, in plasmas 
highly charged ions in a relativistic regime [HID , NEEC has not been observed experimentally yet. 
The theory of NEEC in highly charged ions involves accurate calculations of the atomic structure 
together with a phenomenological description of the nucleus. The free electron and the nucleus 
interact via the Coulomb field (mostly responsible when describing electric multipole transitions 
of the nucleus) or via a current-current interaction (necessary to picture magnetic multipole nu- 
clear transitions). Both the electron-nucleus interaction and the phenomenological description of 
the nuclear dynamics via collective models for NEEC are very similar to nuclear polarization cal- 
culations [ | ll 1|| . The cross section for the two-step process of NEEC followed by the 7 decay of 
the nucleus as a function of the recombined electron energy E can be separated as 

2vr 2 AJT n 



a(E) 



L n (E — E ) 



(28) 



v r n 

where Y n is the NEEC rate and v4 7 the 7 decay rate of the excited nucleus. Furthermore, p is 
the continuum electron momentum, T n the width of the nuclear excited state and L n (E — E ) a 
Lorenzian profile 



L n {E — E ) 



TJ27T 



(29) 



(E - E y + \Yl 

centered on the resonance energy E and with the width given by the width of the nuclear res- 
onance. Due to the very narrow natural width of the nuclear excites states T n , on the order of 
10~ 5 to 10~ 8 eV, the Lorenzian is almost a 5 peak, and the NEEC cross section is non-vanishing 
only very close to the resonance. Since an electron energy resolution on the sub-meV scale is not 
available in present experiments, one should convolute the cross section in Eq. (T2~8T) with a realistic 
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electron energy distribution to simulate experimental data. 

A major difficulty in observing NEEC experimentally arises from the background of atomic 
photorecombination processes, in particular from the non-resonant channel of RR. Due to identical 
initial and final states of NEEC followed by 7 decay of the nucleus and RR, and to the dominance 



of the latter, it is practically impossible to distinguish between the two processes H112H . The 
RR photon yields for typical experimental conditions exceed the ones of the 7 photons following 
NEEC by orders of magnitude, resulting in a signal-to-background ratio of less than 1CT 3 even for 



the most promising cases Illl2h . The convoluted signal to background ratio 



R{E, s) = — — (30) 
a RR (E,s) 

for the Ml nuclear transition from the ground state to the first excited state at £'„=125.358 keV 
with the capture of a free electron into the K shell of a bare 3| 5 Re ion as a function of the free 
electron energy can be seen in Figure [T3J The energy distribution of the incoming electrons is 
assumed to be described by a Gaussian function with the width s. The RR cross section has a 
practically constant value on the energy interval of s. Figure [13] shows the signal to background 
ratio for three different experimental width parameters s = 0.5 eV, 1 eV and 10 eV. While for a 
width parameter s = 0.5 eV the contributions of the NEEC and interference terms can be clearly 
discerned from the RR background, for presently more realistic widths on the order of eVs or tens 
of eV the values of the ratio R(E, s) are too small to be observed experimentally. 

There are at present two ideas how to tackle the problem of observing experimentally NEEC, 
or, to be more specific, how to separate the NEEC signal from the RR background. One of them 
relies on a specific storage ring experimental setup. For NEEC to occur, atomic vacancies are 
needed, i.e., ions. Usually, the higher the ion charge, the larger is the NEEC cross section, such 
that experiments with highly charged ions are desirable. Highly charged ions can be produced in 
accelerator facilities and then transferred and stored in storage rings, where they cycle at almost 
the speed of light. For recombination experiments, the ion beam is driven through an electron 
target. In practice, the different time scales on which NEEC and RR occur can be useful for 
eliminating the background in a storage ring experiment where the ions are cycling in the ring at 
nearly speed of light. While RR is practically instantaneous and its photon signals last typically for 
10~ 14 s, NEEC involves intermediate excited nuclear states which decay after 10~ 10 s or even later 
following the recombination. While the RR photons will be emitted instantaneously in the region 
of the electron target, the radiative decay of the nucleus will occur later, after the ions have already 
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traveled a certain distance in the storage ring. Considering typical experimental parameters of 
the GSI ESR, with the ion velocity being 71% of the speed of light, typical nuclear excited state 
lifetimes of 10-100 ns correspond to 2-20 m distance between the prompt and delayed photon 
emission. The electron cooler is used as a target for free electrons and can be tuned to match 
the resonance condition for NEEC. Ions recombined in the cooler are separated from the primary 
beam in the first bending magnet of the storage ring and can be counted with an efficiency close 
to unity. By choosing isotopes with convenient excited state lifetimes on the order of 10-100 ns, 
the spatial separation can be determined such that a direct observation of the 7 photons following 
NEEC could be performed with almost complete suppression of the RR background 111 1311 . 

Another approach aims at separating the signal and background photons in energy. Typically 
NEEC experiments and theoretical studies envisaged excitation of the nucleus from its ground 
state to the first excited state. However, one can also consider NEEC from a long-lived excited 
state, for instance an isomeric state. Metastable nuclear states (nuclear isomers) are an interesting 
topic in nuclear physics, since it is often not clear why they exi st, and wh ich isomer excitation 
and decay mechanisms are prevailing in different environments ||114l I115Q . In an advantageous 
configuration of the nuclear excited levels (see Figure O, if NEEC occurs from the isomeric state 
to an upper triggering level, the subsequent nuclear decay is possible also to a state below the 
isomer. Such a process is called isomer triggering or depletion, since it allows the depopulation 
of the isomeric state and thus a controlled release of the energy stored in the metastable state. The 
advantage when considering isomer triggering via NEEC is that the energies of the excitation from 
the isomeric state (equal to the RR background photons) and of the photons emitted in the decay 
of the triggering level (which are the signature that NEEC occurred) are different, as shown in 
Figure [141 and thus the experimental background due to RR is much reduced. 

The search for practical methods to trigger isomeric states has been the subject of a number 
of sometimes controversial investigations in the last decades. Major motivation for this was the 
challenge of understanding the formation and the role of nuclear isomers in the creation of the 
elements in the universe. However, the interest has been not least generated by a number of fas- 
cinating potential applications related to the controlled release of nuclear energy on demand, such 
as nuclear batteries that would operate without fission or fusion. A comparison between different 
triggering mechanisms, i.e., different ways to reach the triggering level from the metastable state, 
has surprisingly showed that NEEC is the most efficient way of releasing the isomer energy on 
demand when low-lying triggering levels are used J 1 16 1. These are preferable for obtaining high 
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energy gain and facilitating the excitation to the triggering level. An experimental verification of 
NEEC and isomer triggering via NEEC at the borderline of atomic and nuclear physics would 
be therefore very desirable and will be hopefully rendered possible by upcoming ion storage ring 
facilities and ion beam traps commencing operation in the near future. 



V. CONCLUSIONS 



In atoms, nucleus and electrons are bound together and get to step on each other's toes. The 
electrons are sensitive to the nuclear mass, charge distribution, magnetization distribution, spin, 
and finally also to the nuclear excited levels that can be (virtually) populated. In the search to 
better understand the underlying physics, the last decades have witnessed a substantial increase 
of precision in atomic physics experiments, supported by more and more accurate theoretical cal- 
culations of electronic energies. A summary of the discussed nuclear properties and their effects 
are given in Table U It is usually a challenging task to extract nuclear parameters out of atomic 
physics experimental findings, since most of the time several effects play a role, and each con- 
tribution has to be identified accordingly. Successful approaches often involve putting together 
theory and data obtained with different complementary experimental techniques, many of which 
have been outlined here. 

Strong efforts towards improving both experimental precision and theoretical accuracy in 
atomic and nuclear physics are needed. With the commissioning of new accelerator facilities such 
as the Facility for Antiproton and Ion Research (FAIR) [27] at the GSI, Germany, and the Heavy 



Ion Research Facility in Lanzhou (HIRFL) 111 1711 . China, as well as the increasing number of EBIT 
facilities around the world, many new challenging high-precision experiments with highly charged 
ions are in sight. A very exciting perspective is opened by the Radioactive Ion Beam facilities, 
such as the one of the Michigan State University, USA, TRIUMF in Canada or in RIKEN, Japan, 
which offer the possibility to investigate the atomic and nuclear properties of nuclides far from the 
valley of stability in ionic states. Nuclei not yet available for experiments could be produced in 
quantities that allow determination of their properties. Halo nuclei have already raised a question 
mark on whether our nuclear structure theory at present is applicable only to the limited number 
of nuclei along the stability line. Just as in the case of halo nuclei, atomic physics experiments 
with unstable nuclei might bring new insights in the structure of the nucleus. 

Challenging experiments such as the observation of NEEC or the determination of the nuclear 
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lifetime and properties of the 229m Th isomer will probably benefit of the increasing high precision 
in atomic physics experiments. The 229m Th isomer lies very close to the ground state, the latest 
determination of the excitation energy being 7.6 eV 111 1811 . For the decay of the isomeric state, IC is 
orders of magnitude more probable than 7 decay, and one cannot hope to determine the transition 
energy or the isomer lifetime by detecting the low-energy photon. Instead, one can probe and 
identify the nuclear state by determining the HFS of the electronic levels, taking advantage of the 
different spins and magnetic moments of the ground and isomeric states. The determination of 
the nuclear excited state HFS is demanding but promising for this case, with a number of groups 
around the world engaged in this project. 

Finally, high-precision atomic physics experiments with neutral atoms in parity-violation mea- 
surements on different atoms aim at probing the electroweak interaction and testing the Standard 
Model in search of new physics. A key ingredient is the reliable and accurate theoretical input. 
Whether the atomic theory errors can be reduced to the level of the experimental uncertainty re- 
mains a challenge. An experimental strategy for circumventing the atomic theory uncertainties 
is to measure parity-violating observables for different atoms along an isotopic chain. Standard 
Model predictions for ratios of such observables are largely independent of atomic theory. The 
next decades may show how productive the studies at the borderline between atomic and nuclear 
physics still are and what can atomic transitions reveal about the nucleus, its structure and its 
properties. 
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TABLE I: Nuclear effects in atomic transitions 



Nuclear property 



Effect on atomic structure 



Charge +Ze 

Size - radius trms 

Mass M < oo, nuclear recoil 



Binding energy 
Field shift 

Mass shift: NMS, SMS 



Spin and magnetic moment 
Quadrupole moment 



Magnetic HFS 
Quadrupole HFS 



Weak interaction Parity mixing of atomic states 

Polarizability - virtual nuclear excitations Nuclear polarization 
Nuclear transitions IC/NEEC, BIC/NEET 




12 3 4 5 6 



9 10 



FIG. 1: A two-parameter Fermi nuclear charge distribution determined by the surface thickness t and the 
half-density radius c. 
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600 




X (nm) 



FIG. 2: A typical spectral line obtained from the Is 2s 2p Pi/ 2 — -P3/2 transition in B-like 4U Ar . The 
six dashed curves represent a fit to the Z eem an components. The obtained experimental wavelength for the 
Ml transition is A = 441.2559(1) nm Reprinted with permission from R. Soria Orts et al. , Phys. 

Rev. Lett. 97, 103002, 2006. Copyright (2006) by the Americal Physical Society. 
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FIG. 3: The 7r^ 2 Zeeman components of the ls 2 2s 2 2p 2 Pi/ 2 — 2 Pz/2 transition in B-like 40 Ar 13+ obtained 
from laser spectroscopy at the Heidelberg EBIT. 




1( Li = 9 Li + 2n 



FIG. 4: Artist's view of the halo 11 Li nucleus. The two additional neutrons form the halo around the 9 Li 
core. 
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-fc L — i 8 

-e k — c 



Initial state (i) Resonant state (d) Final state (f) 



FIG. 5: Dielectronic capture of a free electron into the L shell of a He-like atom with the simultaneous 
excitation of one of the core isT-shell electrons. The resonant recombination is followed by the radiative 
decay of the excited electron to the K shell. Such a DR process is called KLL-DR, denoting the involved 
atomic shells. 
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FIG. 6: Experimental data for RR and KL1L2-DR in mercury ions with charge states from He-like to Be- 
like. The emitted photon energy is plotted as a function of the energy of the recombining electron. On top 
of the RR lines one can observe the DR resonances. The principal and total angular momentum quantum 
numbers n and j of the atomic shell into which the electron recombines radiatively are given on each RR 
line. The DR resonances are indexed by j\ and j'2, the total angular momentum quantum numbers of the 
excitation level for the bound electron and the recombination state of the free electron, respectively. 
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Relative Energy [ eV ] 



FIG. 7: Dielectronic recombination of the Li-like neodymium isotopes Nd (orange or gray line) and 
150 Nd 5 (black line) in the energy range of the ls 2 2p 1 / 2 18lj> resonance groups. The f labels indicate the 
individual fine structure components of the n = 18 Rydberg electron. The DR recombination scheme is 
presented in the inset. Reprinted with permission from C. Brandau et at, Phys. Rev. Lett. 100, 073201, 
2008. Copyright (2008) by the American Physical Society. 
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FIG. 8: Illustration of the hyperfme splitting scheme for / = 3/2 and J = 1 considering the example of 
201 Hg. The original fine structure atomic level (left) is split into a magnetic HFS (middle) corrected by the 
quadrupole HFS (right). Compare with Eq. ( fT6l ). 
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FIG. 9: Experimental determination of the F = 4 to F = 3 hyperfine transition in the ground state of 
hydrogenlike 165 Ho 66+ . The signal is compared to the background radiation in the trap when no holmium 
ions are stored. 




FIG. 10: Electroweak interaction between the electron and the nucleus, mediated by the neutral gauge 
boson Z°. 




FIG. 1 1 : The ladder and cross diagrams contributing to nuclear polarization in the lowest order. Here e 
denotes the electron and N and N* the ground state and excited nucleus, respectively. 
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FIG. 12: The inverse process of internal conversion: NEEC. If the atomic and nuclear transition energies 
match, an electron can recombine into an ion (left panel) with the simultaneous excitation of the nucleus 
(right panel) from the ground state G to the excited state E. 
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FIG. 13: The ratio R(E, s) in Eq. (l30l for recombination into bare rhenium as a function of the energy of 



lfl!2 



the continuum electron for three different experimental electron energy width parameters s I112H . See text 
for further explanations. 
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FIG. 14: Partial level scheme of ||Mo. The isomeric state (/) can be excited (for instance via NEEC) to the 
triggering level (T) which subsequently decays back to / or to a level F, initiating a cascade via different 



intermediate states (dashed line) to the ground state (GS) 1112011 . The direct I —t F decay is a strongly 



hindered EA transition, while / — > T and T — > F are E2 transitions. 
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